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Abstract: The level-crossing probability, local and global adiabaticity conditions are discussed for 
2-flavour neutrino oscillations in matter with arbitrary mixing angles "d. Different approximations 
for the survival probability of supernova neutrinos are compared. Results of a combined likelihood 
analysis of the observed SN 1987A neutrino signal and of the latest solar neutrino data including the 
recent SNO CC measurement are presented. 



1. Neutrino evolution: resonance 
and adiabaticity conditions, max- 
imal violation of adiabaticity 

We consider neutrino oscillations in a two flavour 
scenario and label the heavier neutrino mass 
eigenstate with "2" . Then A = m\ — m\ is pos- 
itive and the vacuum mixing angle i? is in the 
range [0 : 7r/2]. As starting point for our dis- 
cussion, we use the evolution equation for the 
medium states vb 




where A m = {(A - Acos2z9) 
denotes the difference between the effective 
masses of the two (active) neutrino states in 
matter, E is their energy and A = 2EV = 
2^/2GpN e E is the induced mass squared for the 
electron neutrino. Furthermore, $ m is the mix- 
ing angle in matter and <&' m — di9 m /dr. Since 
anti-neutrinos feel a potential V with the oppo- 
site sign than neutrinos, formulae derived below 
for neutrinos become valid for anti-neutrinos af- 
ter the substitution i9 — > n/2 — 

The traditional condition for an adiabatic 
evolution of a neutrino state along a certain tra- 
jectory is that the diagonal entries of the Hamil- 
tonian in Eq. are large with respect to the 
non-diagonal ones, |A m | ^> \4E"3' m \. This con- 
dition measures indeed how strong adiabaticity 
is locally violated. Therefore, the point of max- 



imal violation of adiabaticity (PMVA) is given 
by the minimum of A m /?9( n . In the following, 
we will concentrate on power-law like potential 
profiles, A cx r n . This type of profile does not 
only contain the case n ss — 3 typical for super- 
nova envelopes, but also the exponential profile 
of the sun in the limit n — > ±oo. Moreover, it 
allows discussing which features of neutrino os- 
cillations are generic and which ones are specific 
for the linear profile n = 1 usually discussed. For 
A(r) oc r™, the minimum of A m /d / m is at 



cot(2tf, 
-- [cot(2tf 



2tf) +2cot(2i? m ) (1.2) 
2<d) - cot(2tf m )] = . 



For n — 1, the PMVA is indeed at i? m = 7r/4 for 
all ■&. Thus, in the region where the resonance 
point flm = 7r/4 is well-defined, they coincide. 
In the general case, n ^ 1, the PMVA agrees 
however only for •& — with the resonance point. 

In Fig. 0, we show the the change of the sur- 
vival probability, Ap(r)/dr = dl^^P/dr for a 
neutrino produced at r — as 9%, together with 



the PMVA predicted by Eq. (1.2) and the res- 



onance point for a power law profile A oc r . 
The resonance condition predicts a transition in 
lower-density regions than the PMVA, until for 
= 7r/4 the resonance point reaches r = oo and 
the concept of a resonant transition breaks down 
completely. Moreover, the crossing probability 
becomes less and less localized near the PMVA 
for larger mixing angles 
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Figure 1: Change of the survival probability 
dp(r)/dr of a neutrino produced at r = as V2 to- 
gether with the point of maximal violation of adia- 
baticity (dot) and the resonance point (star) for a 
power law profile A oc r~ 3 . The height of the differ- 
ent curves is rescaled. 

Let us now discuss the condition for the adi- 
abatic evolution of a neutrino state along a tra- 
jectory from the core of a star to the vacuum. 
While the condition |A m | > |4P$JJ indicates 
whether adiabaticity is locally violated, we need 
now a global criterion that measures the cumu- 
lative non-adiabatic effects along the trajectory 
from -d nl ir/2 to For a non-adiabatic evolu- 
tion of the neutrino state we require that 









/ dtf m Vi 


= e 


/ dtf m 


Jtt/2 




J IT j2 



(1.3) 



with e< 1. Then the border between the adia- 
batic and non-adiabatic regions is given by 



2n{2V Q ) 



A _ _ 

E ~ \ £ iTn 2 (2i?)(l - sintf) Ra 
with 



l/n ^ n+1 



(1.4) 



tt/2 



dtf m sintf m [sin(27? m )] 2+1/ " 
x [sin(2 l ? m -2t9)] 1 - 1/n |. (1.5) 



Fig. H shows the excellent agreement between 
our prediction for the border between the non- 
adiabatic and adiabatic regions for anti-neutrinos 




W 10 s 
tan 2 (i?) 



Figure 2: Contours of constant anti-neutrino sur- 
vival probability (dashed) together with the border- 
line Eq. between adiabatic and non-adiabatic 
regions using Eq. (1.5) with n = —3 (solid) and 



n — > ±oo (dash-dotted) in /(#) for the SN profile 
given in the text; the dotted line shows the border- 
line for neutrinos. 

with energy E = 20 MeV and a profile typical 
for supernova envelopes, V(r) = 1.5 x 10~ 9 eV 
(10 9 cm/r) 3 , and the one following from the con- 
tours of constant survival probability P ee (dashed 
lines) of the neutrino eigenstate v-i obtained by 
solving the Schrodinger equation (1.1). A com- 
parison of the solid (n = —3) and the dash-dotted 
line (n — ► ±oo) shows moreover that /(??) de- 
pends only weakly on n. 

2. The crossing probability in the 
WKB formalism 

The leading term to the crossing probability 
Plsz within the WKB formalism is in the ultra- 
relativistic limit and omitting an overall phase 
given by 



In Plsz = - 



1 

E 



dx A, 



x 1 (A 1 ) 



(2.1) 



where A2 = iAe 2 " 5 are the branch points of A m 
in the complex x plane and and X2 can be chosen 
arbitrarily either on the positive or negative real 
x axis. The usual choice, A\ = AC, allows to 
express In Plsz as the product of the adiabatic- 
ity parameter 7 evaluated at the resonance point 
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and a correction function T n (2) that can be rep- 
resented as a hypergeometric function 2 Pi [[|. 

Another representation for the crossing prob- 
ability which is valid for all ■& uses as integra- 
tion path in the complex x plane the part of a 
circle of radius A centred at zero and starting 
from Ai = A and ending at A 2 = Ae 2li} . For 
A = A (r I Ro) n , one can factor out the 1? depen- 
dence of Plsz into functions Q n 0, 



where 



]nPj jSZ = -K n GnW, (2-2) 



is independent of and 



9? / d^e^ (e inv -C) +S 



(2.4) 




with C = cos2tf and 5" = sin2tf. 



Figure 3: Comparison of the contours of constant 
survival probability Pee calculated numerically for a 
M = 20A/q progenitor star (dotted lines) and cal- 
culated for A oc r -3 with the Plsz approximation 
(solid lines). 



3. Neutrino oscillations in super- 
nova envelopes 

In the analysis of neutrino oscillations, the poten- 
tial profile A(r) of supernova (SN) envelopes is 
often approximated by a power law with n « —3, 
and V(r) = 1.5 x 1(T 9 eV (10 9 cm/r) 3 . A com- 
parison of the results of a numerical solution of 



the Schrodingcr equation (1.1) with the analyt- 
ical calculation of P gg using the G-3 functions 
shows very good agreement for this profile; only 
tiny deviations in the region A/(2E) ~ 1CT 17 eV 
have been found in 0. By contrast, all other ap- 
proximations used hitherto in the literature fail 
in some part of the tan 2 i?-A plane: while the use 
of T\ = 1 together with A cx r~ 3 describes cor- 
rectly the crossing probability for small mixing 
in the resonant region, the errors become larger 
for larger mixing until this approximation fails 
completely in the non-resonant region. The cor- 
rection function Too used for n = — 3 describes 
quite accurately the most interesting region of 
large mixing as well as the non-resonant region, 
but does not reproduce the correct shape of the 
MSW triangle. 

More important is however to check how 
strong deviations of the true SN progenitor pro- 
file V(r) from a power-law profile may affect 



the analytical results. Realistic progenitor pro- 
files differ in two aspects from a simple 1/r 3 be- 
haviour. First, the outer part of the envelope has 
an onion like structure, and its chemical compo- 
sition, Y e (r), and thus also V(r) changes rather 
sharply at the boundaries of the various shells. 
Second, the density drops faster in the outermost 
part of the envelope, becoming closer to an expo- 
nential decrease. We calculated numerically Pgg 
using profiles for different progenitor masses and 
stellar evolution models. We found that Pgg is 
well approximated only in the non-resonant part 
by our analytical results for the 1/r 3 profile, inde- 
pendently of the details of the progenitor profile, 
while Pgg depends strongly on the details of the 
progenitor profile in the resonant region. As an 
example, we compare in Fig. ^the Pgg calculated 
numerically for a 20M© profile with the analyt- 
ical results for our standard SN profile. There- 
fore a numerical solution of the Schrodinger equa- 



tion (1.1) should be performed in the resonant 
region, using a realistic profile for the particu- 
lar progenitor star considered. However, a 1/r 3 
profile together with the WKB crossing probabil- 
ity is sufficient for the analysis of tmfci-neutrino 
oscillations in the phenomenologically most in- 
teresting region tan 2 •& < 5 independent of the 
details of the progenitor envelope. 
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Next we present the results of a combined 
statistical analysis of the neutrino signal of 
SN 1987A and of the complete set of solar neu- 
trino experiments 8. Since the two data sets 
are statistically independent and functions of the 
same two fit parameters, they can be trivially 
combined, 

Xt 2 ot (tf, A) = X 2 (tf, A) + x§ N (0, A) . (3.1) 

Contours of constant confidence level (C.L.) are 
defined relative to the minimum of Xtoti where 
Xq(i9, A) was calculated in Ref. || for the solar 
data and X sn = -2£(i?,A) in Ref. @ for the 
SN 1987A data. We consider the astrophysical 
parameters as known and minimize only the two 
parameters ■& and A. 

In Fig. [I] we show the C.L. contours of the 
combined fit for a rather representative set of as- 
trophysical parameters, namely binding energy 
E h = 3 x 10 53 erg and (E Qe ) = 14 MeV. In this 
case, the impact of the SN 1987A data on the 
standard solutions to the solar neutrino prob- 
lem is rather dramatic: the LOW-QVAC and 
VAC solutions disappear for both assumed r = 
(Ep h )/(E Pc ) values; they are excluded at more 
than 99.98% even for r = 1.4. Moreover the size 
of the LMA-MSW solution decreases with in- 
creasing r. The part of the LMA-MSW solution 
which is most stable against the addition of the 
supernova data corresponds to the lowest A and 
tan 2 i? values, since these are favoured by Earth 
matter regeneration effects. On the other hand 
the SMA-MSW region re-appears extending, for 
increasing r, as a funnel towards the VAC solu- 
tion along the hypotenuse of the solar MSW tri- 
angle. The combined best-fit point (star) moves 
from the LMA-MSW region for r = 1.4 to the 
SMA-MSW solution for r = 1.7. 

In Fig. [| we illustrate in a global way the 
relative status of various oscillation after adding 
the SN 1987A data. For each A value we have 
optimized the \ 2 with respect to ?9 in the top 
and with respect to A in the bottom panel. The 
solid curve indicates Xqj the non-solid curves 
correspond to the case where the SN 1987A 
data are included. The dash-dotted line is for 
E h = 3 x 10 53 erg, r = 1.4 and (E 0<s ) = 14 MeV. 
The dashed line is for E h = 3 x 10 53 erg, r = 1.4 
and (Eg e ) = 12 MeV. The dotted line is for E h = 




10 I 

Iff" 10' 3 Iff 2 10'' 1 10 

tan 2 (tf) 
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Figure 4: The 90, 95, 99 and 99.73% C.L. con- 
tours of the combined fit of solar and SN 1987A data 
(coloured/grey) together with the contours of the so- 
lar data alone (solid lines); for r = {Ea h ) / (Ea e ) = 1.4 
(top) and r = 1.7 (bottom). All figures for Eh = 
3 x 10 53 erg and {Es> B ) = 14 MeV. 



3 x 10 53 erg, r = 1.7 and (E Pe ) = 14 MeV. Here 
we have adjusted an arbitrary constant which ap- 
pears when combining the minimum likelihood- 
type SN 1987A analysis with the solar x 2 data 
analysis in such a way that the SMA solution gets 
unaffected by the SN 1987A data. One notices 
that the effect of adding SN 1987A data is always 
to worsen the status of the large-mixing angle so- 
lutions. Within each such curve one can compare 
the relative goodness of various solutions, how- 
ever different curves should not be qualitatively 
compared. 



International Europhysics Conference on HEP 



M. KachelricB 



80 



70 



60 



50 - 



40 - 



30 




10""10" 10 10-" 10' 8 10~ 7 10-" 10" 5 10' 4 10" 3 

A(eV 2 ) 




Figure 5: The solid curve indicates the x© °f 
the solar neutrino data. The non-solid curves illus- 
trate the effect of adding the SN 1987A data, which 
worsens the status of large mixing-type solutions; 
marginalized with respect to ■& (top) and to A (bot- 
tom), respectively. 



solar neutrino data including the SNO CC mea- 
surement, we found that the supernova data of- 
fer additional discrimination power between the 
different solutions of the solar neutrino puz- 
zle. Unless all relevant supernova parameters 
lie close to their extreme values found in sim- 
ulations, the status of the LMA solutions dete- 
riorates, although the LMA-MSW solution may 
still survive as the best combined fit for accept- 
able choices of astrophysical parameters, fn par- 
ticular, SN 1987A data generally favour its part 
with smaller values of $ and A. In contrast 
the vacuum or "just-so" solution is excluded and 
the LOW solution is significantly disfavoured for 
most reasonable choices of astrophysics parame- 
ters. The SMA-MSW solution is absent at about 
the 3<7-level if solar data only are included but 
may reappear once SN 1987A data are added, 
due to the worsening of the LMA type solutions. 

Finally, one should not forget that in the so- 
lar well-tested standard solar model ex- 
ists whose errors are accounted for in the fit. In 
contrast there is no "standard model" for type II 
supernovae and therefore also no well-established 
average values and error estimates for the rele- 
vant astrophysical parameters. 
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4. Summary and Discussion 



We have discussed non-adiabatic neutrino oscil- 



lations in general power-law potentials A oc x" . 
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